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ABSTRACT 

We model an accretion disk atmosphere and corona photoionized by a central X-ray continuum source. 
We calculate the opacity and one-dimensional radiation transfer for an array of disk radii, to obtain the 
two-dimensional structure of the disk and its X-ray recombination emission. The atmospheric structure 
is extremely insensitive to the viscosity a. We find a feedback mechanism between the disk structure 
and the central illumination, which expands the disk and increases the solid angle subtended by the 
atmosphere. We apply the model to the disk of a neutron star X-ray binary. The model is in agreement 
with the ~ 12° disk half-angle measured from optical light curves. We map the temperature, density, 
and ionization structure of the disk, and we simulate high resolution spectra expected from the Chandra 
and XMM-Newton grating spectrometers. X-ray emission lines from the disk atmosphere arc detectable, 
especially for high- inclination binary systems. The grating observations of two classes of X-ray binary 
systems already reveal important spectral similarities with our models. The model spectrum is dominated 
by double-peaked lines of H-like and He- like ions, plus weak Fe L. The line flux is proportional to the 
luminosity and is dominated by the outer radii. Species with a broad range of ionization levels coexist 
at each radius: from Fe XXVI in the hot corona, to C VI at the base of the atmosphere. The line 
spectrum is very sensitive to the temperature, ionization, and emission measure of each atmospheric 
layer, and it probes the heating mechanisms in the disk. We assume a hydrostatic disk dominated by 
gas pressure, in thermal balance, and in ionization equilibrium. As boundary conditions, we take a 
Compton-temperature corona and an underlying Shakura-Sunyaev disk. The choice of thermally stable 
solutions strongly affects the spectrum, since a thermal instability is present in the regime where X-ray 
recombination emission is most intense. 

Subject headings: accretion, accretion disks — atomic processes — instabilities — line: formation — 
X-rays: binaries 



1. INTRODUCTION 

When the infall of matter into a deep gravitational po- 
tential is mediated by an accretion disk, gravitational en- 
ergy is converted to the thermal radiation which powers 
both low mass X-ray binary (LMXB) systems and ac- 
tive galactic nuclei (AGN). Accretion disks present unique 
problems involving magnetized plasma dynamics, pho- 
toionization, atomic kinetics, thermal and ionization equi- 
libria, general relativity, and radiation transfer. The ac- 
cretion disks in LMXBs and AGN are expected to have 
many common properties. The compactness of the ac- 
cretor in LMXBs and AGN and their inferred accretion 
rates imply temperatures of T > 10 7 K and intense X- 
ray emission in the inner disk region. The inner radii of 
these disks, as well as the disk atmosphere as a whole, is 
substantially more ionized than the case in which the ac- 
cretor is a white dwarf. In both LMXBs and AGN, the 
vast energy emitted in the inner disk region is reprocessed 
in the outer disk, where the external radiative heating can 



dominate the local thermal emission. The subsequent pho- 
toionization of the disk plasma radically alters its equilib- 
rium state, structure, and spectrum, especially in the at- 
mospheric and coronal disk layers, which are the subject 
of this study. 

High resolution X-ray spectroscopy is an essential tool 
to study the physics of this "hot class" of accretion disks 
and the conditions near black hole event horizons. In this 
paper, we concentrate on the outer radii of disks in neu- 
tron star LMXBs, since current observational constraints 
provide more stringent tests for LMXBs than for AGN. 
The following points support these assertions: 



High resolution spectra can reveal discrete emission 
or absorption from atomic transitions within 
the accretion disk plasma, providing information 
on the accretion disk structure, dynamics, and 
physics. These spectra open a window into 
photoionized gases and their phase equilibria. 
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• X-rays, and in particular discrete atomic tran- 
sitions of hydrogen- and helium-like ions, probe 
the regions in the disk with the highest levels of 
ionization. Regions closer to the compact object 
will have the highest ionization levels, although 
vertical stratification is also expected. 

• The knowledge of the accretion disk physics 
directly impacts our ability to probe the physical 
conditions around the compact object. For 
example, the Fe K emission originating in the 
innermost regions of an accretion disk has 
been proposed as a direct probe of the general 
relativistic effects near a black hole event horizon 
in AGN, by virtue of the observed characteristic 
line shape (Tanaka et al. 1995). However, very 
little is known about the physical conditions in the 
Fe K emission region, and it is still unclear how 
our ignorance of the physical processes within the 
accretion disk affect the modeled Fe K line profile 
and flux. It is also unclear whether the soft X-ray 
line features reported by Branduardi-Raymont et 
al. (2001) are feasible. 

• While in neutron star LMXBs the photoionizing 
source must be near the neutron star surface, in 
AGN and galactic black hole candidates (BHC) 
the location of the ionizing source is unknown. In 
AGN, various authors have assumed the ionizing 
source to be located in the rotation axis of the 
black hole, above the disk midplane, possibly close 
to the base of a jet. Alternatively, an ionizing 
source might be present on the upper layers 

of the disk, perhaps due to disk flares, or to 
Comptonization of thermal UV photons in the 
accretion disk corona (ADC). 

• LMXB systems are observed in less crowded 
regions than AGN. 

• In contrast to AGN, LMXBs often have measured 
orbital parameters which constrain the geometry 
of the system, such as the maximum disk radius. 
LMXBs may also have a measured value of the 
disk inclination, while orbital phase and eclipse 
phase variations provide tomographic information. 
For example, measurements of the optical light 
curve amplitudes of LMXBs have yielded estimates 
of the angle subtended by the disk of ~ 12° (de 
Jong, van Paradijs, & Augusteijn 1996). 

• To our knowledge, the inner disk radii can only 
be studied in the X-ray band. Theoretically, the 
thermal emission of the inner disk radii must peak 
in the X-ray band for both LMXBs and AGN. The 
optical and UV emission originates in the outer 
regions of the accretion disk and further away from 
the compact object, and the radio and gamma ray 
emission is likely dominated by emission from jets. 

Fewer physical ingredients are needed to model the outer 
radii of disks than the inner radii, so the logical progres- 
sion is to successfully model the outer radii first. The X- 
ray spectroscopy of neutron star LMXBs will allow us to 



construct a physical picture of their accretion disks, which 
we can then use to investigate the inner disk in BHC and 
AGN. 

To fully exploit the high energy-resolution X-ray spec- 
tra of accretion disks, we created physical disk models and 
calculated synthetic spectra. The disk plasma, at 10 5 to 
10 7 K, cools through atomic line emission that can be de- 
tected with space-borne X-ray observatories such as Chan- 
dra and XMM-Newton. Modeling the equilibrium state 
of the plasma and the radiation transfer within the disk 
allows a calculation of the disk structure and its X-ray 
spectrum. A synthetic spectrum can be compared to the 
data. The model spectrum is unique in that it is calcu- 
lated purely on physical, and not just phenomcnological, 
grounds. 

We describe four fiducial disk models. Two of these 
models were introduced in Jimenez-Garate, Raymond, 
Liedahl, & Hailcy (2001). We use a newly developed 
adaptive-mesh disk structure calculation, the Raymond 
(1993) photoionized plasma code, and a new X-ray emis- 
sion code which uses HULLAC data (Hebrew Univer- 
sity/Lawrence Livermore Atomic Code, Klapisch et al. 
(1977)). The models consist of a disk illuminated by a 
pure neutron star continuum, and they contain as bound- 
ary conditions a Compton-temperature ADC at the top of 
the disk and a modified Shakura & Sunyaev (1973, here- 
after abbreviated as SS73) disk at the bottom (Vrtilek et 
al. 1990). Thus, the region of interest has a temperature 
and ionization which is intermediate of these two regions, 
and it emits copious X-ray radiation. 

In section 2.1, we introduce the X-ray line observations 
prior to Chandra and XMM-Newton; in section 2.2, we 
introduce theoretical work on the structure of X-ray illu- 
minated accretion disks; in section 3, we describe the disk 
structure calculations and the assumptions of hydrostatic, 
thermal and ionization equilibrium; in section 4, we detail 
the effects of a thermal instability on a layer of the disk 
atmosphere throughout the disk; in section 5, we discuss 
the calculation of the high resolution spectrum, which is 
done a posteriori from the structure calculation; in section 
6, the disk density, temperature and ionization structure 
are presented; in section 7, the model spectra are shown, 
assuming a full, partial or obstructed view of the neu- 
tron star region, and we show simulated spectra utilizing 
the response of the XMM-Newton reflection grating spec- 
trometer (RGS) and the Chandra medium energy gratings 
(MEG); in section 8, comparisons of the model to the ob- 
served X-ray spectra of LMXBs are discussed briefly, and 
we discuss the limitations of the model. In section 9, con- 
cluding remarks are presented. 



2. LMXB ACCRETION DISKS 

2.1. X-ray line emission from LMXBs 

With the exception of Fe K emission in the 6.4-7.0 keV 
range (Asai, Dotani, Nagase, & Mitsuda 2000), discern- 
ing X-ray line emission in LMXBs has been challeng- 
ing, owing to limitations in sensitivity and spectral re- 
solving power, as well as the difficulties associated with 
attempts to extract line emission from data dominated 
by intense continuum emission. Measurements obtained 
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with the Einstein Objective Grating Spectrometer (Vr- 
tilek et al. 1991), the ROSAT Position Sensitive Propor- 
tional Counter (Schulz 1999), and the ASCA CCD imag- 
ing detectors (Asai, Dotani, Nagase, & Mitsuda 2000) 
have shown that the spectra of a large fraction of bright 
LMXBs exhibit line emission. X-ray lines at ~ 1 keV are 
often mixed with various species, so that only the bright- 
est LMXBs had clear line identifications, as in the case of 
Ne X Lya in 4U1626-67 (Angelini et al. 1995), Fe L in Sco 
X-l, or O VIII Lya & Ly/3 in 4U1636-53 (Vrtilek et al. 
1991). 

The X-ray line emission arises presumably as the result 
of irradiation of the disk by the X-ray continuum, pro- 
ducing an extended source of reprocessed emission. Evi- 
dence of X-ray emission from extended regions in LMXBs 
comes from the spectral variations during ingress and 
egress phases of eclipses, and during rapid intensity fluc- 
tuations known as dips. Most dips, which are observed to 
precede eclipses, are thought to result from variable ob- 
scuration and attenuation of the primary continuum by 
material near the outer disk edge, which has been thick- 
ened due to impact of the accretion stream with the disk 
(White & Swank 1982; Frank, King, & Lasota 1987). Dips 
that are uncorrelated with orbital phase can be produced 
by orbiting clouds crossing the line of sight, as shown in 
Figure 1. A cloud larger than ~ 10 6 cm can obscure the 
X-rays from the neutron star. Hard X-ray emission, pre- 
sumably originating in the ADC, and representing a few 
percent of the non-eclipse flux, remains visible during mid- 
eclipse in several LMXBs, implying that the ADC is larger 
than the secondary star (White & Holt 1982; McClintock, 
London, Bond & Grauer 1982). LMXB spectra during 
eclipses or dips may harden or soften, i.e. the propor- 
tion between hard (<~ 3-10 keV) and soft (~ 1-3 keV) X- 
rays changes. Most sources harden during dips (Parmar, 
White, Giommi & Gottwald 1986), consistent with photo- 
electric absorption, but there are exceptions like the soft- 
ening of 4U1624-49 (Church & Balucinska-Church 1995), 
and an unchanging X1755-33 (White et al. 1984; Church 
& Balucinska-Church 1993). Sources such as EXO0748- 
67, X1916-05, and X1254-69 show evidence for an unab- 
sorbed spectral component during dips (Parmar, White, 
Giommi & Gottwald 1986; Church et al. 1997), reveal- 
ing an extended source of X-rays which is larger than the 
ADC. These soft X-rays are likely radiation reprocessed 
in the accretion disk. Dip ingress/egress times indicate 
ADC sizes in the 10 9 to 5 x 10 10 cm range, a factor of a 
few smaller than the accretion disk sizes calculated from 
typical orbital parameters (Church 2001). A soft X-ray 
emission component distinct from the hard X-ray contin- 
uum has also been interpreted as being due to the effects of 
absorption edges or line emission in some LMXBs (Parmar 
et al. 2000; Church et al. 1998). The unequivocal identifi- 
cation of the accretion disk X-ray emission, requires both 
an energy resolution that is higher than CCD detectors, 
and high throughput, plus a quantitative theoretical pre- 
diction of the X-ray emission from the disk. We discuss the 
recently observed high resolution spectra in section 8.1. 



2.2. Radiatively heated accretion disks 
2.2 A. Radial structure 



In LMXBs roughly half of the gravitational potential en- 
ergy is released in the vicinity of the compact object (i.e., 
in a boundary layer near the neutron star surface). The 
disk is exposed to this radiation, and it will be heated by 
it. Radiative heating can exceed internal viscous heating 
in the outer region of the disk. The temperature structure 
of the disk can thus be controlled by the X-ray field pho- 
toionizing the gas, suppressing convection, and increasing 
the scale height of the disk. 

Assuming that all the viscous heating and radiative 
heating from illumination by the central source is radiated 
locally as a blackbody (as in the SS73 model), Vrtilek et 
al. (1990) calculated the temperature for a geometrically 
thin disk, with r 3> Ri, where R\ is the radius of the 
compact X-ray source: 

SGMiM + (l-r))L x sm0(r) 



phot ~ 87rr 3 ' 4irr 2 
where T p hot is the photospheric temperature, M\ is the 
mass of the compact X-ray source, G is the gravitational 
constant, a the Stcphan-Boltzmann constant, 9 is the graz- 
ing angle of the incident X-ray flux with respect to the 
disk surface, and r\ is the X-ray albedo such that (1 — 77) 
is the fraction of X-rays absorbed at the photosphere. 
The albedo has been derived from optical observations (de 
Jong, van Paradijs, & Augusteijn 1996). The first term on 
the right-hand side of equation (1) is the energy dissipated 
within the SS73 disk, and the second term is the radiative 
heating. The radiative heating term will dominate where 



r > 2.3x10* 



Mi 



1 



V 



0.1 



sinb 

~oT 



0.1 



cm, 
(2) 

where the X-ray luminosity is written in terms of an X- 
ray accretion efficiency e x , according to L K — e x Mc 2 . For 
example, accretion onto a neutron star results in roughly 
1/2 of the gravitational potential energy being converted 
into X-rays, or e x = GM\f2c 2 R\. The disk, therefore, is 
radially divided in two regions: an inner region dominated 
by internal dissipation, and an outer region dominated by 
external illumination. External radiation will dominate 
the disk atmosphere energetics for the outer two or three 
decades in radii, and the local dissipation and magnetic 
flare heating, if any, will be ignored there (see also section 
8.3). 

2.2.2. Vertical Structure 

The radial dependence of the disk temperature in equa- 
tion (1) relies on averaging physical quantities such as the 
dissipation parameter a in the direction perpendicular to 
the disk plane, which is valid for regions in the disk that 
are optically thick. However, as we will show in this article, 
the radiative recombination spectrum is very sensitive to 
the radial and the vertical ionization structure, including 
regions with an optical depth r < 1. 

To obtain a high resolution spectrum of an accretion 
disk, and in particular one for which the outer (or upper) 
layers are X-ray photoionized, several authors have calcu- 
lated the vertical structure by solving the radiation trans- 
fer equations, assuming hydrostatic equilibrium. Models 
have been applied to AGN and LMXBs in the high-L x 
state, since in the low-state radiatively inefficient accretion 
ensues, which is described by a separate family of models 
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(Hawley & Balbus (2002) and references therein). In ra- 
diativcly efficient accretion disks, the radiative transfer is 
typically simplified by using an on-the spot approxima- 
tion and the escape probability formalism. Due to photo- 
electric absorption and Compton scattering, the ionization 
structure of the disk is stratified, and it is approximated 
by a set of zones, each with a single ionization parameter. 
The ionization structure of the disk can be solved by us- 
ing photoionization codes such as CLOUDY (Ferland et 
al. 1998) and XSTAR (Kallman & McCray 1982), which 
calculate the ionization and thermal equilibrium state of 
the gas at each zone. 

Ko & Kallman (1991, 1994) calculated the vertical struc- 
ture of an illuminated accretion disk and obtained the re- 
combination X-ray spectrum for individual rings on the 
disk. Raymond (1993) utilized the temperatures in equa- 
tion (1) and calculated the vertical structure and the UV 
spectrum from the entire disk. Both assumed parame- 
ters for LMXBs, and gas pressure-dominated disks. Later 
models of photoionized accretion disks focused primarily 
on calculating the Fe Ka fluorescence emission from AGN 
disks. 

Rozahska & Czerny (1996) and Rozahska, Czerny, Zycki 
& Pojmahski (1999) modeled semi-analytically the strat- 
ified, photoionized transition region between the corona 
and the disk in AGN. They found that their approxi- 
mations, which included on-the-spot absorption, matched 
more accurate radiation transfer codes for optical depths 
< 10. They also discussed the existence of a two-phase 
medium, stopping short, however, of calculating an X-ray 
spectrum. Nayakshin, Kazanas & Kallman (2000) mod- 
eled a radiation-pressure dominated disk and showed that 
the vertical structure of the disk implied significant dif- 
ferences in the Fe K fluorescence line spectrum compared 
to that predicted by constant-density disk models (Ross 
& Fabian 1993; Matt, Fabian & Ross 1993; Zycki, Kro- 
lik, Zdziarski & Kallman 1994). In addition, Nayakshin, 
Kazanas & Kallman (2000) also found that the gas was 
thermally unstable at certain ionization parameters, which 
created an ambiguity in choosing solutions and a sharp 
transition in temperature in the disk. This instability is 
discussed in section 4. Ballantyne, Ross, & Fabian (2001) 
calculated the vertical structure of disk ring as a func- 
tion of radius, accretion rate, the angle of incidence of 
radiation, the photon index, and the black hole mass, al- 
beit using a diffusion approximation. Rozahska, Dumont, 
Czerny, & Collin (2002) calculated the hydrostatic disk 
structure including Compton scattering and line transfer 
without assuming the escape probability approximation. 
Rozahska, Dumont, Czerny, & Collin (2002) also calcu- 
lated the structure of the optically thick part of the disk, 
by use of the diffusion approximation and the local a- 
prescription (Eq. [12]). All of the above models calculate 
the disk structure for one radius at a time. 

Li, Gu, & Kahn (2001) found the static solution that 
resolves the thermal instability in the gas by considering 
the effect of conduction, and they computed the X-ray 
recombination and resonance line scattering spectrum for 
the conduction transition region that forms between stable 
solutions in the disk. With this procedure, the unphysical, 
sharp transition between stable phases was eliminated. Li 
et al. considered ionizing continua typical of AGN, which 



can yield stable solutions with three different tempera- 
tures for a given pressure ionization parameter S (defined 
in section 3.3). Up to three distinct transition layers can 
form. The reflection and recombination spectrum of the 
transition regions in the 0.5-1.5 keV range was computed 
by considering the vertical structure of an isobaric, opti- 
cally thin region. They found that resonant scattering can 
be important within the transition region, depending on 
the local gravity and luminosity, which yields a line spec- 
trum which is different from that of pure recombination 
emission. 

The vertical structure of an optically thick accretion 
disk can be obtained using the diffusion approximation, 
which assumes that the photon mean free path A is much 
smaller than the scale of temperature and density gra- 
dients T/VT and p/Vp, respectively. Adding convective 
heat transfer by introducing an adiabatic temperature gra- 
dient, Meyer & Meyer-Hofmeister (1982) have calculated 
the vertical structure of an isolated accretion disk which 
is dominated by convection zones. Such techniques are 
used in the standard stellar structure equations. X-ray 
illumination from the central compact object suppresses 
convection, reduces the thermal gradients in the disk, and 
has a stabilization effect in the outer radii; but if X-ray il- 
lumination is combined with the diffusion approximation, 
it also produces a convex disk that self-shadows the outer 
disk regions. This contradicts the observed spectra, which 
show evidence of reprocessing from the outer disk (Dubus 
et al. 1999). A semi-analytical model using a variable a- 
viscosity prescription was used to model AGN disks and to 
investigate its effects on the Lyman edge absorption and 
emission (Rozahska, Czerny, Zycki & Pojmanski 1999). 

The failure of the diffusion-equation models to repro- 
duce a concave disk that can efficiently reprocess the cen- 
tral X-rays may indicate that important effects were ne- 
glected. First, the effect of the disk atmosphere and corona 
was ignored. Second, turbulent heat transfer may pro- 
duce a vertical disk structure that is nearly isothermal. 
The strong turbulence occurring at the scale of the disk 
thickness in magncto-hydrodynamic (MHD) models sup- 
ports this hypothesis (Miller & Stone 2000). A reliable 
calculation of the turbulent heat transfer in an accretion 
disk is needed. Therefore, we prefer to use the Vrtilek et 
al. (1990) vertically-isothermal disk for the optically thick 
region. 

The diffusion approximation is inadequate when calcu- 
lating high resolution spectra since line radiation must 
originate in a region where the photon mean-free-path A 
exceeds the scale of the temperature gradient, i.e. A > 
T/VT. Thus, just as for stellar atmospheres (Mihalas 
1978), an explicit radiation transfer calculation without 
assumption of local thermodynamic equilibrium is needed. 
The modeling of photoionization heating, recombination 
cooling, and X-ray opacities is then required in the atmo- 
sphere. 

3. MODEL ATMOSPHERE 

We consider a LMXB with a M* = 1.4 M & primary ra- 
diating an Eddington luminosity (L x = 10 38 - 3 erg s _1 ) 
bremsstrahlung continuum, with T — 8 keV. A set of 
fiducial system parameters for a bright LMXB is used, 
so application to a particular source will require using the 
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observed X-ray continuum to improve accuracy. The max- 
imum radius of the centrally-illuminated disk is 10 11 cm, 
so the orbital period ~ 1 day. The minimum radius is 
10 8 ' 5 cm, below which the omitted effect of radiation pres- 
sure, in large part, determines the disk structure. 

The vertical structure of the disk atmosphere for each 
annulus in the array, is obtained by integrating the hydro- 
static balance and 1-D radiation transfer equations for a 
slab geometry (Fig. 2): 

dP GM*pz 

la w 



dz 
dF v 

dz 

dFl 

dz 



sm( 



(4) 



(5) 



while satisfying local thermal equilibrium (see also eq. 
[16]): 

A(P,p,F v ) = (6) 

and ionization balance (see also eq. [13]): 

ion formation rate = ion destruction rate (7) 

where P is the total pressure, p is the mass density, F v is 
the net flux of incident radiation (which is the intensity 
integrated over all solid angles in erg cm~ 2 s _1 Hz -1 ), 
is the reprocessed net flux propagating down towards the 
disk midplanc, z the vertical distance from the midplane, 
G the gravitational constant, 6 the grazing angle of the ra- 
diation on the disk, r the radius, v is the frequency, and n v 
is the local absorption coefficient. The rays corresponding 
to F v and F* are defined in Figure 3. Hydrostatic equilib- 
rium is satisfied to a < 1% accuracy and thermal balance 
to < 0.01%. 

For the structure calculation only, 100 logarithmically 
spaced energy bins, in the range 1 eV < hv < 1000 keV, 
were used for F v and F^. The grid is coarse, and yet suf- 
ficiently broad to accomodate a hard X-ray tail in future 
models. The reprocessed radiation propagating upwards, 
F", is omitted to accelerate the computation. This is a 
good approximation since the radiative heating is domi- 
nated by the direct flux F v . The reprocessed flux F" is 
calculated a posteriori by a high resolution spectral model 
(section 5). The difference between cooling and heat- 
ing, A, includes Compton scattering, bremsstrahlung cool- 
ing, photoionization heating, collisional line cooling, and 
recombination line cooling (section 3.4). Cosmic abun- 
dances (Allen 1973) are assumed. The code from Ray- 
mond (1993) computes the net heating and ionization 
equilibrium, models Compton scattering in one dimen- 
sion, and calculates line scattering using escape proba- 
bilities. A new disk structure calculation simultaneously 
integrates equations (3)-(5) by the Runge-Kutta method, 
using an adaptive step-size control routine with error es- 
timation, and equation (6) is solved by a globally conver- 
gent Newton's method (Press 1994). At the ADC height 
z COI , the equilibrium T is close to the Compton temper- 
ature T compt on, from which we begin to integrate down- 
ward until T < TphotM- The optically thick part of the 
disk, with temperature Tphot, is assumed to be vertically 
isothermal (Vrtilek et al. 1990). To get Tphot, the viscous 
energy and the illumination energy are assumed to be lo- 
cally (re)radiated with a blackbody spectrum. Thus, for 



Zphot ^ r an d R* ^ r , equation (1) can be used with 
Mi = M*. The height at which the integration ends is 
defined as the photosphere height z p hot- Thus, we assume 
that viscous dissipation dominates heating for z < z p hot 
(Fig. 3). 

The boundary conditions, shown schematically in Fig- 
ure 3, are set at the ADC to P(z COI ) = p C orkT coinpton / p,m p , 
J F v (z cor )dv — L x /4irr 2 , and F^(z cor ) = 0, where k is the 
Boltzmann constant, and p is the average atomic weight 
of baryons in units of the proton mass m p . The bound- 
ary conditions at the photospheric height (= Zphot) for 
F v and F^ are set free, and the shooting method (Press 
1994) is used with shooting parameter p cor , which is ad- 
justed until P(z p hot) — Pphot&Tphot/ pm p is satisfied at the 
photosphere. Note p p hot is the viscosity-dependent density 
calculated for an X-ray illuminated SS73 disk. 

The shooting method consists of guessing the value of 
the coronal density which matches the desired pressure at 
the bottom of the gas column. The boundary conditions 
define P cor once p COI is chosen. Equations (3)-(7) are si- 
multaneously solved during the integration. The tempera- 
ture drops as the integration proceeds downward through 
the atmosphere. When T reaches a value below Tphot, the 
pressure at that point is compared to the expected pres- 
sure of the isothermal disk at that height. If it does not 
match to better than ~1 %, the integration is repeated 
with a new estimate of the coronal density. While it is not 
clear that photoionization will cease to be important for 
temperatures less than Tphot, such zones emit negligible 
X-ray fluxes if r > 10 9 cm. Our new structure calculation 
also includes the effects of physical instabilities (section 
4), and it removes the numerical instabilities obtained by 
Raymond (1993). 

A novel and important feature of this model is that the 
incident radiation is allowed to modify the disk atmosphere 
geometry, such that the heating and expansion of the at- 
mosphere resulting from illumination are used to calculate 
the height profile of the atmosphere as a function of ra- 
dius. This feedback between the radiative heating and the 
atmospheric structure is depicted in Figure 4. The atmo- 
spheric height is used to derive the input grazing angle 
of the radiation for the next model iteration. This con- 
trasts with calculating the grazing angle using the pressure 
scale height of the optically thick disk (Vrtilek ct al. 1990), 
which is in general much smaller than the photoionized at- 
mosphere, and which underestimates the grazing angle and 
the line intensities by an order of magnitude. To get Tphot 
self-consistently from equation (1), the equation 



0(r) ~ (3 — a + arctan — 



= arctan ^ ^ — arctan^ 



^atn 



+ 



aii lan( — 

r 



(8) 



is needed, where z atm (r) is defined as the height where 
the frequency-integrated grazing flux J F^(z at m) dp is at- 
tenuated by e _1 , and u,f3 are defined in Figure 2. The 
arctan(F*/r) term is neglected, which is valid for r > 
10 8 ' 5 cm. As discussed above, 8(r) is calculated iteratively 
from equation (8). After an initial guess for z atm (r), it is 
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re-calculated from the newly obtained disk structure. A 
power-law fit to z atm (r) works well to obtain 9(r). This it- 
eration is performed with a limited number of radial bins 
(5), to save computation time. The iteration is stopped 
after 8(r) and T p hot converge to < 10 %. After conver- 
gence, the number of logarithmically spaced radial bins is 
increased to 26. The process of convergence does not de- 
pend on the initial choice of 0(r), and it is shown in Fig- 
ure 5. However, convergence does depend on the choice of 
Zatmi which is a free parameter in the model. Since z a tm is 
not physically determined, it must be defined ad hoc, but 
it is bound by z atm > z pho t- For z < 

z phot, the illumina- 
tion F v — > 0, and the disk blackbody flux Fbb (Tphot) takes 
over. To test how sensitive is the result to the definition 
of Zatm, we calculate 8{r) taking z atm = z p hot, and we use 
this to estimate the systematic errors of the 1-D radiation 
transfer calculation. 

3.1. The Choice of Assumptions 

The validity of the assumptions is reviewed, both for the 
the model presented here and for some of the accretion disk 
models in the astrophysical literature. 

For modeling X-ray line emission from the disk atmo- 
sphere, the commonly used assumptions of LTE and the 
diffusion approximation will not hold. In addition, assum- 
ing a constant density in the vertical direction will be inad- 
equate, since the hydrostatic equilibration time is small or 
comparable to other relevant timescales, and the line emis- 
sion is highly sensitive to the vertical ionization structure. 
The recombination emission is especially sensitive to this 
structure, since each ion emits clearly resolvable line en- 
ergies. Also, fluorescence emission can be reprocessed by 
a Compton-thick, fully ionized gas above it (Nayakshin, 
Kazanas & Kallman 2000). 

Thermal equilibrium and ionization equilibrium are rea- 
sonable assumptions in a time averaged sense. Hydrostatic 
equilibrium is also assumed, to avoid explicit computation 
of the plasma dynamics with radiative transfer. Devia- 
tions from hydrostatic equilibrium are smoothed in the 
timescale 

^hydro = Z atm /C s ~ 2 - 7 ^ 1 Q7 t c m ) 7 5~ 1/2 SCC ( 9 ) 

where c s is the sound speed. Material from the disk 
moves radially within the viscous timescale i v i sc ~ 
r 2 /(az atm c s ) > thydro (Frank, King & Raine 1992), so the 
gas can reach hydrostatic equilibrium before it flows in- 
ward (i.e, the radial accretion velocity is always subsonic). 
However, the Kcplerian orbital velocity 3> c s is highly 
supersonic. If the gas flow in the corotating frame of the 
gas is also supersonic, then shocks would collisionally ion- 
ize and heat the gas. In such a case, the observed spectrum 
of the disk would significantly deviate from a photoionized 
gas in ionization and thermal equilibrium. The thermal- 
ization and ionization equilibrium timescale of the atmo- 
sphere is driven by the recombination timescale 

/ j,l/2 s 

tth = t rec ~0.3 — sec (10) 

which can be derived from equation (A3), and where T5 
is the temperature in units of 10 5 K, nu is the den- 
sity in units of 10 14 cm~ 3 , and Z is the atomic number 



(Reynolds & Fabian 1995). The photoionization timescale 
is shorter than i roc where the gas is fully stripped; other- 
wise both timescales are comparable for the relevant ions 
having similar abundances. The Coulomb collision relax- 
ation timescale between electrons and ions t cp is slower 
than between identical particles, and is (Spitzer 1962) 

rp3/2 

t C p ~ 3 x 10~ 6 ^— sec. (11) 

ni4 

Electron-electron relaxation is ~ m p /m c times faster, and 
proton-proton relaxation is ~ \J m p /m c times faster. For 
the hot corona at the outer disk at T ~ 10 7 K and 
n c ~ 10 11 cm~ 3 (from the coronal structure in section 6), 
the relaxation timescale is t cp ~ 3 sec. The fully ionized 
gas in the corona, which is near the Compton temperature, 
has a thermal timescale of tth — ^compton = lOr^o^is sec > 
where no is the disk radius in units of 10 10 cm, and £38 
is the X-ray luminosity in units of 10 38 ergs" 1 (Reynolds 
& Fabian 1995). Thus, thcrmalization in the disk atmo- 
sphere and corona is driven by the ionization timescales, 
since the Coulomb relaxation times are comparatively fast 
due to the large density. Thermal and ionization equilib- 
rium occur faster than hydrostatic equilibrium, for length 
scales 2: atm > 10 6 cm. If t t h < ihydro, luminosity fluctua- 
tions with timescales ta ux such that tth < iflux < ihydro will 
take the gas outside hydrostatic equilibrium, but not out 
of thermal equilibrium. Integrated spectral observations 
on timescales t 3> tfl U x cannot observe this effect. 

The radiation transfer is complex, and the assumptions 
used to simplify calculations could be problematic. In par- 
ticular, by dividing the disk atmosphere into annular zones 
with a given vertical gas column, our 1-D radiation trans- 
fer calculation assumes that 1) the primary continuum is 
not absorbed before reaching the top of the column, 2) the 
radiation in the column propagates from top to bottom at 
a given grazing angle, and 3) there is no significant radia- 
tive coupling from one disk annulus to another, which is 
used to justify the slab approximation. The above assump- 
tions are inadequate if the column height is comparable to 
the disk radius, or if the photon mean free path in the gas 
column is many times the local radius. Thus, future 2-D 
calculations will result in better bookkeeping of photons, 
a more accurate structure, and a more reliable X-ray line 
spectrum. 

The correct calculation of line transfer in the gas is also a 
concern, since the disk atmosphere is optically thick in the 
lines (section 6). Line transfer is complicated by the Keple- 
rian velocity shear, which has to be taken into account for a 
given viewing angle (Murray & Chiang 1997). The escape- 
probability approximation used to calculate line transfer 
in the disk may also be inadequate because of the large 
optical depths. 

Our calculations show that the proper treatment of a 
thermal instability (Field 1965; Krolik, McKee & Tarter 
1981) and conduction affect the spectrum significantly 
(Zeldovich & Pikelner 1969; Li, Gu, & Kahn 2001) (sec- 
tion 4). A two-phase gas could form, with clouds of an 
unknown size distribution and with undetermined dynam- 
ics of evaporation and condensation (Begelman & McKee 
1990), with each phase having distinct ionization param- 
eter and opacity. The instability is sensitive to 1) the 
metal abundances, 2) the continuum shape (Hess, Kahn, 
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& Paerels 1997), and 3) the atomic kinetics (Savin et al. 

1999) . 

The local viscous energy dissipation rate per unit vol- 
ume in the disk atmosphere may be included in equation 
(6) with the form (Czcrny & King 1989, SS73) 

Q visc = ^ftaP (12) 

where f2 is the Keplerian angular velocity a is the viscos- 
ity parameter, and P is the local pressure. Equation (12) 
is an extension of the a-disk model (where the viscous dis- 
sipation is vertically averaged), and it assumes the local 
validity of the a prescription, which is untested. Fortu- 
nately, our numerical modeling indicates that the viscos- 
ity term is negligible in most regions of the disk atmo- 
sphere except for the inner disk and a ~ 1 (in particu- 
lar, near the Compton-temperature corona). This viscos- 
ity term enhances a thermal instability between 10 6 and 
10 7 K. Vertically stratified MHD models (Miller & Stone 

2000) , although inconclusive, owing to the uncertain effect 
of boundary conditions, show that the viscous dissipation 
drops rapidly at > 2 pressure scale heights away from the 
disk midplane, providing evidence against equation (12). 
Since the disk atmosphere is always a few scale heights 
above the midplane, we choose not to include equation 
(12) in our models. Equation (12) has been applied in the 
optically thick regions of the disk by Dubus et al. (1999) 
and in the disk atmosphere by Rozahska & Czerny (1996). 
Other forms for the local dissipation that reduce to the 
a-disk have been used (Meyer & Meyer-Hofmeister 1982). 

3.2. Ionization balance 

In steady state the equation of ionization balance for 
each ion Z +l is 

= n z ^ i+1 n c a Zti+1 + n z . i ^ 1 (f3 z ^ 1 + n e C Z)i _i) 

-n z ,i{(3 Z! i + n c a Ztl + n e C z ,i) 
+n z , l - 2 [3l l _ 2 B z . l _ 1 = 0, (13) 

where (3 z .i is the photoionization rate (s _1 ) of Z +l , C z .i is 
the collisional ionization rate coefficient (cm 3 s" 1 ) of Z +l , 
a Zi i + \ is the recombination rate coefficient (cm 3 s" 1 ) of 
ion Z + ( l+1 \ and so on. The code also includes K-shcll 
photoionization followed by Auger ionization in the last 
term of equation (13). The Auger branching fraction is 
B z ,i-i = 1 — Yz,i~i, where Y z ^-\ is the fluorescence yield. 
Multiple Auger decays are ignored in equation (13). The 
terms with a z ^ + i and ot z ^ account for all two-body recom- 
bination processes. The coefficients C z ,i and a z .i depend 
on the electron temperature T for any ion Z +t . Given 
the photoionization cross-section o~pE,z,i and ionization 
threshold energy \z,i of the ion Z +t , the photoionization 
rate for a point source of ionizing continuum is 

= £ j^^^W^M (14) 

where Se is the spectral shape function, normalized on 
a suitable energy interval. For the accretion disk atmo- 
spheres orbiting neutron stars, which are of interest here, 
collisional ionization rates are negligible compared to pho- 
toionization rates. 



3.3. The ionization parameter 

Let £ = L x /n p r 2 (in erg s" 1 cm) be the ionization pa- 
rameter, where n p is the proton number density (Tarter, 
Tucker & Salpeter 1969). The ionization parameter £ is 
factored out of equation (14), and together with the spec- 
tral shape function Se, it defines uniquely the charge state 
distribution in an optically thin photoionized gas (eq. [13] 
does not include three-body recombination, which is im- 
portant at n e > 10 16 cm -3 ). 

Another, dimensionless ionization parameter is con- 
structed with the radiation pressure P ra d and the proton 
gas pressure P gas (Krolik, McKee & Tarter 1981). This 
ionization parameter, S, is defined as S e Prad/ Pgas, 
where P rad = J F v dv/c. Note S = ^/inckT. The new 
parameter S is useful when the local pressure can be de- 
fined. For an optically thin gas, an isobar has constant 

3.4. Thermal Equilibrium 

We review the terms in the thermal equilibrium equa- 
tion, which we solve with the Raymond (1993) photoion- 
ized plasma code. Thermal equilibrium is enforced at each 
zone in the disk atmosphere. The explicit form of the ther- 
mal equilibrium condition, equation (6), is 

Compton net heating + photoionization heating = 

bremsstrahlung cooling + recombination cooling 

+collisional cooling, (15) 

which corresponds to (Halpern & Grindlay 1980) 

«e5> M (A£r + K Ti + a c 4) (i6) 

in units of erg s" 1 cm" 3 , where the rates for each process 
and other dependencies are included in the group of rate 
coefficients (r™ m , T^, A h z rcm , A^ c , A c °j), n is the elec- 
tron number density, n z ^ is the Z +l ion density, and the 
sums are performed over all abundant ions. The radiative 
heating is directly proportional the net flux F v (in units 
of erg s" 1 cm" 2 Hz" 1 ), and the density. Cooling pro- 
cesses, which originate from electron-ion interactions, are 
proportional to the square of the density and are in general 
dependent on the electron temperature T. The coefficients 
in equation (16) can be obtained from Halpern & Grindlay 
(1980), and some, such as the recombination coefficients, 
are very dependent on the available atomic data. A list 
of the data used for the coefficients in the model, and a 
list of the processes and transitions included in the calcu- 
lations can be found in Raymond (1993). Recombination 
cooling includes both radiative recombination and dielec- 
tronic recombination. Collisional cooling includes cooling 
due to line emission and collisional ionization. Photon 
trapping and subsequent collisional de-excitation reduces 
the cooling rate compared to the optically thin case, but 
this affects the UV lines more than the X-ray lines because 
the resonant scattering opacity is larger in the UV. 

The ions of H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, and 
Fe are included in the thermal equilibrium equation (16), 
and the ionization balance equation (13). 
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For a fully ionized gas, such as the hot corona above 
the accretion disk, Compton heating and inverse-Compton 
cooling dominate equation (16). The Compton net heat- 
ing term may be positive or negative, since the transfer of 
energy between the photons and the electron gas depends 
on the electron temperature and the shape of the ionizing 
spectrum. In such cases, and in the non-relativistic case, 
the equilibrium temperature is 

h J vF v dv 
Tcompton = 4kjF u dv 

where k and h are the Boltzmann and the Planck con- 
stants, respectively. The Compton temperature T compton 
is determined uniquely by the shape of the ionizing contin- 
uum. For an 8 keV bremsstrahlung spectrum, T compt on ~ 
2 x 10 7 K. 

4. THERMAL INSTABILITY IN PHOTOIONIZED GASES 

Irradiated gas is subject to thermal instabilities for tem- 
peratures in the 10 5 -10 7 K range (Buff & McCray 1974; 
Field 1969), such that X-ray line emission at those tem- 
peratures may be suppressed. The Field (1965) stability 
criterion, together with plasma equilibrium calculations 
(Davidson & Netzer 1979; Kallman & McCray 1982; Ray- 
mond 1993; Ferland et al. 1998), indicate that a photoion- 
ized gas becomes thermally unstable when recombination 
cooling of H- and He-like ions is important. Consequently, 
the disk atmosphere structure has a thermally unstable 
region, which modifies the X-ray spectrum. 

To clarify the nature of the thermal instability, consider 
the calculated net heating (Fig. 6). The gas is externally 
heated, and the net heating depends on the state variables 
and ionization of the gas, making this a peculiar system. 
The thermal balance locus, where the net heating is zero, 
is denoted as the S-curve and is displayed in Figure 6. The 
region to the left of the S-curve undergoes net cooling, and 
the one to the right has net heating. 

To test stability, consider small T-perturbations start- 
ing from the S-curve. A vertical displacement from any 
point in Figure 6 represents an isobaric perturbation in 
T. From this, we find that the points in the S-curve with 
positive slope are stable, while those with negative slope 
are unstable (Field 1965). This splits the S-curve into 
three branches. The shape of the S-curve determines which 
range of S are unstable (section 3.3). At such S, thermal 
balance is achieved by three distinct T on the S-curve, two 
stable T and one unstable T. On the unstable branch, iso- 
baric T-perturbations cause a thermal runaway to one of 
the stable branches. 

The shape of the S-curve depends on the metal abun- 
dances and the ionizing spectrum (Hess, Kahn, & Paerels 
1997). The S-curve is subject to uncertainties in the 
atomic data (Savin et al. 1999), and its shape may vary 
(albeit not dramatically) from one plasma code to another. 
Our calculated S-curve for the disk atmosphere is shown 
in Figure 7. 

Most spectral studies of heated accretion disks in 
LMXBs and AGN have either used unstable solutions, or 
just selected a subset of the stable solutions. Rozahska, 
Dumont, Czerny, & Collin (2002) chose a monotonic den- 
sity, which is equivalent to selecting all points on the S- 
curve. This results in a pressure which oscillates with 



height and a transition region which is not in hydrostatic 
equilibrium. Ko & Kallman (1994) and Nayakshin & Kall- 
man (2001) selected the hot branch of solutions, which 
produce a condensing atmosphere biased towards high- 
ionization species. Ballantyne, Ross, & Fabian (2001) do 
not specify how the choice of solutions within the insta- 
bility was made, but they acknowledge the effects of the 
instability, which are seen in the sharp temperature tran- 
sition obtained with their models. 

The instability implies a large VT as the gas is forced to 
move between stable branches, requiring the formation of a 
transition region whose size may be determined by electron 
heat conduction, convection, or turbulence, depending on 
which dominates the heat transfer. For simplicity, calcu- 
lations of emission from the transition region are omitted 
in this article. Upon calculation of the Field length \p, 
the lengthscale below which conduction dominates thermal 
equilibrium (Begelman & McKee 1990), we estimate that 
conduction forms a transition layer ~ 10~ 2 times thinner 
than the size of the X-ray emitting zones. Nevertheless, 
X-ray line emission from the neglected conduction region 
may not be negligible in all cases (Li, Gu, & Kahn 2001). 
The £ values present in the transition region are absent 
in other regions, which may allow the transition region to 
have observable spectral signatures. Resonant scattering 
from the transition region may be observable in some sit- 
uations (Li, Gu, & Kahn 2001). The importance of the 
transition region also depends on the shape of the S-curve 
and on the local gravity. 

Conduction tips the balance of stability at sufficiently 
small spatial scales. If the gas is not in static equilibrium, 
conduction can drive phase transitions, and produce dy- 
namic condensing or evaporating fronts. In static equilib- 
rium, conduction quenches the instability and produces a 
transition layer at S sta t (stretching the S-curve in Fig. 7) . 
This transition layer (Zeldovich & Pikelner 1969) connects 
the low-T stable branch at S < S sta t with the high-T sta- 
ble branch at 5 > S sta t- In the dynamic case, if the tran- 
sition layer is located away from S stat , it will dynamically 
approach S sta t by a conduction driven mass flow, as shown 
in Figure 8. A transition layer with S sta t < H < S evap pro- 
duces an evaporating front, while a transition layer with 
S CO nd < S < S sta t produces a condensing front (Zeldovich 
& Pikelner 1969; Li, Gu, & Kahn 2001). 

The disk structure for both condensing and evaporating 
solutions is computed. We assume a steady state, condens- 
ing or evaporating mass flow through the transition layer 
at S con( j or S evap , respectively. The static equilibrium so- 
lution is an intermediate case of the latter extreme cases. 
A single-valued T(S) is used, since a two-phase solution 
may be buoyantly unstable, making the denser (colder) 
gas sink. The evaporating disk corresponds to the low-T 
branch, while the condensing disk corresponds to the high- 
T branch (Figure 7). This introduces spectral differences 
(section 7). 

We do not know from first principles whether the disk 
atmosphere is evaporating, condensing, or static. How- 
ever, a Compton-heated wind might be expected in the 
corona for large radii (Begelman, McKee & Shields 1983). 
The speed of the conduction mass flow is estimated to be 
Wcond = 2KT/3P gas AF, by using the characteristic conduc- 
tion time at the Field length, where k is the Spitzer (1962) 
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conductivity (for a detailed discussion, see McKee & Begel- 
man 1990). A conduction mass flow speed w con d = l - 
2xl0 -2 times the local sound speed is obtained. Thus, the 
phase dynamics will depend on the subsonic (v > w CO nd) 
flow patterns in the disk atmosphere, and these flows will 
in part determine the evaporation or condensation rates, 
together with the boundary conditions on mass flow. 

If the disk is in a steady state of evaporation or con- 
densation, the implied mass flow can have an effect on 
the global mass budget, due to mass conservation. Steady 
state evaporation implies mass loss or a disk wind, while 
condensation implies a mass gain (Zeldovich & Pikclner 
1969; Li, Gu, & Kahn 2001). 

A thermal instability due to Compton heating and 
bremsstrahlung cooling can ensue between 10 6 and 10 7 K 
if the ionizing spectrum extends well above ~ 10 keV (Kro- 
lik, McKee & Tarter 1981). For a 8 keV bremsstrahlung 
spectrum, this additional instability regime is suppressed 
(Hess, Kahn, & Paerels 1997). Nevertheless, some LMXBs 
have harder spectra (White et al. 1995). A double S-curve 
results from the hard spectra in AGN (Nayakshin & Kali- 
man 2001), which allows a three-phase gas. 

As mentioned above, gas dynamics which are not in- 
cluded in the model can have an impact on the gas phase. 
The only physical mechanism known to transport the nec- 
essary angular momentum for disk accretion involves a 
magneto-rotational instability (MPJ) which drives turbu- 
lent flow in the disk (Balbus & Hawley 1998). These tur- 
bulent flows are nearly supersonic in the disk midplane re- 
gion, where most of the mass is accreted. Enhanced heat 
transfer rates due to this turbulent flow could quench the 
thermal instability and affect the disk structure. Turbu- 
lent heat transfer rates can be orders of magnitude larger 
than the saturated conduction heat transfer rate. How- 
ever, it is not known whether such turbulent motions will 
also be present in the disk atmosphere, which is several 
scale heights above the disk midplane and has a density 
which is orders of magnitude smaller (Section 6). A de- 
cline in the viscous a parameter with vertical disk height 
was obtained with local MHD models, and an enhanced 
ratio of the magnetic pressure to the gas pressure with 
increasing height (Miller & Stone 2000). The MRI also 
favors the assumption of vertical isothermality in the op- 
tically thick disk. In section 8.3, we discuss the effects of 
magnetic fields in the atmosphere and corona. 

5. SPECTRAL MODELING 

With the disk structure p(r,z), T(r,z) and ion abun- 
dances f z .i(r, z), the X-ray line emission from the disk at- 
mosphere is calculated using HU LLAC data (Klapisch et 
al. 1977). The code calculates the atomic structure and 
transition rates of radiative recombination (RR) and the 
ensuing radiative cascade, which can produce both line 
photons and radiative recombination continuum (RRC) 
photons. We include the H-like and He- like ions of C, 
N, O, Ne, Mg, Si, S, Ar, Ca, and Fe, as well as the Fe L 
shell ions. Fluorescence emission, which is prominent for 
high-Z ions such as those of Fe, is omitted in these calcu- 
lations, as well as resonant scattering, an additional source 
of line emission. The recombination emissivities, and the 
opacities in this model, are calculated as described in ap- 
pendices A and B, respectively. 



The spectrum for each of the 26 annuli was added to 
obtain the disk spectrum. Each annulus consists of a grid 
of zones in the vertical z direction, and T, p and f z ,i+i 
for each zone are used to calculate the RR and RRC emis- 
sivities. The radiation is propagated outwards at inclina- 
tion angle i, including the continuum opacity of all zones 
above, thus accounting for the optical depth of the atmo- 
sphere. Compton scattering of the irradiating continuum 
is included in the disk structure calculation but it is omit- 
ted in the synthetic spectrum. The latter scattering adds 
a weak continuum component with the spectral shape of 
the neutron star emission. The spectrum is Doppler broad- 
ened by the projected local Keplerian velocity, assuming 
azimuthal symmetry. 

6. DISK STRUCTURE 

Once the atmosphere and corona are accounted for, the 
disk is thicker than would be expected from the local pres- 
sure scale height (= Zp) alone. We found Zp < z p hot < 
Zatm- To quantify the disk geometry, the calculated height 
of the photosphere and atmosphere, z p hot and z a t m , are 
both fitted with z = C(r/1 cm)™, with fit parameters C 
and n. The fitted parameters are C p hot = (2-4^°;g) x 10~ 3 
cm, n phot = 1.14±g-g?, C7 atm = (l.Otft?) x 10~ 3 cm, n atm = 
1.21 ±0.01. The above fits imply z phot ~ 3Z P to AZ P (de- 
pending on radius) and z atm ~ (7±l)Z P to (8± 2 )Zp. We 
account for statistical errors and estimated systematics. 
Vrtilek et al. (1990) estimated n atm = 9/7 = 1.29, but in 
spite of the steeper radial dependence, the Vrtilek et al. 
disk is thinner, and it assumes z atm = Zp for r > 10 10 cm. 
The disk thickness derived from the optical light curve ob- 
servations of LMXBs relies on the large fraction of X-rays 
from the neutron star which are shielded from the com- 
panion by the disk. This de facto disk boundary should 
be taken to be ~ z atm , since, by definition, a fraction 1/e 
of the central X-rays are absorbed there. In Figure 9, we 
compare Zp with z a tm- We find that previous theoretical 
studies severely underestimated the size of the disk atmo- 
sphere. 

The X-ray continuum opacity of the atmosphere is 
t C 1 for most lines of sight, except for rays orig- 
inating on the neutron star which are incident at a 
small grazing angle 9(r), such that they are nearly 
parallel to the disk plane. The atmosphere's (z > 
z P hot) maximum photoelectric opacity is always r <C 
1 in the vertical direction, although t/ sin9(r) > 
1. Illumination heating dominates at r > 10 10 cm 
(cq. [1]), where only / F„(z p hot)<W / F v (z COT )dv ~ 
0.12 of the incident photons reach z p hot directly, while 
/ F^(zph.ot)du/(J F u (z COT )dvsm6) ~ 0.35 reach z pho t after 
reprocessing in the atmosphere. Thus, the atmospheric 
albedo is ~ 0.5. Both the photosphere and atmosphere 
contribute significantly to the disk albedo. The total disk 
albedo deduced by de Jong, van Paradijs, & Augustcijn 
(1996) from optical observations in LMXBs is r\ ~ 0.9. We 
have found that its high value is partially explained by the 
atmospheric contribution. Once the latter is taken into 
account, the photosphere's albedo becomes ~ 0.8, since 
0.5+0.5(.8)=0.9, which is closer to physical expectations. 

At any fixed radius, the vertical disk structure has a 
marked boundary between the optically thick, colder disk 
and an optically thin, hotter atmosphere, as shown in Fig- 
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ures 10 and 11. At the largest scales, the vertical struc- 
ture has two distinct zones: a hot corona, in which Comp- 
ton heating and cooling dominates, and an atmosphere or 
warm corona, where photoionization heating and recom- 
bination cooling are most important. Three regions are 
discernible in Figure 12, which shows the vertical struc- 
ture of the outer radius of the disk (other radii show a 
similar pattern, aside from changes in scale). 

The structure of the underlying atmosphere is better 
discerned by plotting the height of the atmosphere above 
the photosphere, z — z p hot- This reveals the presence of 
fine structure, in particular a region emitting lines from 
low-Z He-like ions at T ~ 5 x 10 4 K (Fig. 13 and 14). 
The evaporating and condensing disk model solutions (sec- 
tion 4) are shown in Figures 10 through 14. This low-Z 
He-like ion region is small due to the rapid increase in 
continuum opacity with decreasing temperature, but is 
resolved by the adaptive step-size integration. A more 
extended, T ~ 10 6 K region emits predominantly H-likc 
ion and mid-Z He-like ion RR lines. Both H-like and He- 
like ion emission regions can be identified by the abun- 
dance distribution of the fully ionized and H-like ions, 
which rccombine to produce the H-likc and He-like ion 
emission, respectively (Fig. 15 and 16). The recombina- 
tion line luminosity for a u — > I transition in ion Z +% is 
dL u ^i oc n c n z . l+1 T-~<dV (eq. [A9] and eq. [A10]). The 
highest emissivities will be produced at low temperatures 
and high densities. This implies the region of origin of the 
emission will track the abundances from Figures 15 and 
16, with an added bias towards the lower range of heights, 
which are denser and colder. 

The spatial distribution of K and L-shell Fe ions shows 
that the structure calculation included all the intermedi- 
ate ionization states (Fig. 17). The ionization parameter 
S varied continuously with atmospheric height from full 
ionization at S ~ 10 3 , down to the thermal instability 
regime at S ~ 10, where a break occurs. 

The presence of the instability has a large effect on the 
luminosity of He-like ion lines from mid-Z elements, as 
can be seen from Figure 16. In particular, the Mg +n 
abundance is never allowed to peak, such that the model 
predicts a dim Mg XI line. A similar effect occurs with 
Si XIII and Ne IX. Thus, in the context of this disk model, 
the brightness of these three lines will determine whether 
the instability is operating as modeled. 

The discontinuity in the density, temperature, and ion- 
ization state is a result of enforcing the thermal stabil- 
ity of the chosen solutions, since a range of temperatures 
from ~ 6 x 10 4 K to ~ 7 x 10 5 K is unstable (see section 
4). The discontinuity is unphysical, of course, and can 
be smoothed in future models by the inclusion of conduc- 
tion, or any other heat transport mechanisms in the disk 
that might be present, such as those due to turbulence or 
convection. 

A comparison of the spatial ion distribution from the 
condensing disk (Fig. 18) and the evaporating disk (Fig. 
15) shows that the differences in the synthetic spectra can 
be attributed to differences in the vertical disk structure. 
The low-Z Hc-likc ion line producing region shrinks, while 
the H-like ion line and mid-Z He- like ion line producing re- 
gion expands in the condensing disk model, as compared 
to the evaporating case. The condensing solution shows 



a more extended Fe L emission region, with particularly 
large ion emission measure for Fe +19 , which recombines 
to Fe +ls and produces strong Fe XIX lines (Fig. 17[b]), 
while the evaporating disk shows a larger emission mea- 
sure for Fe +17 at lower temperatures (Fig. 17[c]), which 
recombines to Fe +16 and emits the Fe XVII lines more effi- 
ciently, as will be shown in section 7. This behavior traces 
back to the choice of solutions from the stability curve in 
section 4. 

7. SPECTROSCOPY 

In this section we delineate the circumstances under 
which the disk emission is rendered observable, and we 
describe the disk spectroscopy and its diagnostics. 

The LMXB photon net flux (photon cm" 2 s" 1 keV" 1 ) 
is modeled by: 

F| ot = e- aEN »F* E + e - aEN " ek F E isk (18) 
where F E is the neutron star continuum, F E lsk is the RR 
line and RRC modeled flux (Fig. 19), N£ and N$ ak are 
the neutral hydrogen absorption column densities, E is the 
photon energy, and ge are the Morrison & McCammon 
(1983) absorption cross sections. The system is assumed 
to be d = 10 kpc away. 

We find that the lines are swamped by the continuum for 
cases where the neutral column densities for the neutron 
star and the disk are set equal, or N* H = Nfj sk in eq. [18] 
(see the spectrum in Fig. 20 [a]). This situation is most 
likely to occur in LMXBs with inclination in the range 
i = 0°-60° (Frank, King, & Lasota 1987). The inclination 
angle i is defined in Figure 1. Moreover, the continuum X- 
ray emission from the inner disk (r < 10 8 5 cm) has been 
neglected here, which according to a model by Stella & 
Rosner (1984), will soften the continuum below ~ 10 A, 
and this will further reduce the equivalent widths of the 
X-ray emission lines from the outer disk. In the inner disk, 
radiation pressure dominates and the SS73 viscosity pre- 
scription must be modified (Stella & Rosner 1984). Thus, 
low inclination neutron star LMXBs are unlikely to have 
detectable X-ray lines from the disk. 

Thus, consider N H = 5 x 10 22 cm" 2 and N$ ak = 10 21 
cm~ 2 with i — 75°, where an obscuring medium absorbs 
half of the continuum flux from the neutron star. Such a 
medium is compact enough to leave the disk almost un- 
obscured (see Fig. 20 [b]). This situation should arise in 
LMXBs which exhibit flux dips, for example, where cither 
the disk rim or small clouds obscure the central continuum 
periodically, as explained in section 2.1. These LMXBs 
have inclinations in the 60° to 80° range (Frank, King, 
& Lasota 1987). With a partially obscured central con- 
tinuum, disk evaporation has an observable spectral sig- 
nature. We simulated 50 ks observations with the XMM- 
Newton RGS 1 and the Chandra MEG (Fig. 21 and 22). 
Some bright lines are listed on Table 1. The evaporating 
and condensing disks have contrasting O VII/O VIII and 
Ne IX/Ne X line ratios. The evaporating disk contains 
gas at T ~ (7 TO) x 10 4 K, unlike the condensing disk. 
The H-like ion line intensities are higher for the condens- 
ing disk since it has more gas at T ~ 10 6 K. The spec- 
tral differences stem from the distinct differential emission 
measure distributions d(EM) jd log S and from the O VII 
recombination rate a^R oc T~ 7 , where 7 = 0.7-0.8, and 
EM = J n z ^n e dV is the emission measure (appendix A). 
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In the case the central continuum is completely occulted, 
the model predicts that numerous hard X-ray lines will 
become detectable with Chandra (Fig. 20 [c]), such that 
evaporating and condensing disk models are distinguish- 
able. Figure 23 contrasts the Chandra MEG +1 simula- 
tions of the evaporating and condensing disk models. The 
column density for the neutron star continuum is taken as 
Nfj = 10 24 cm" 2 , and A$ sk = 10 21 cur 2 for the disk. 
Notably, the He-like to H-likc ion line ratios still serve 
to differentiate the models at larger Z, but with the re- 
verse effect. The Hc-like/H-like ion line ratios for Ar, S, 
and Si are larger for the hotter, condensing disk. The 
Mg XI/Mg XII ratio is roughly the same for either model. 
The Q =(Fe XXV + Fe XXVI)/Si XIV line ratio is 50 % 
larger for the evaporating disk model. Since Fe XXV and 
Fe XXVI lines originate near or at the hot corona, where 
the instability in question does not operate, no difference 
in their line fluxes is observed. Since the hot atmosphere 
(or "warm corona") of the condensing disk is larger than 
the evaporating case, more Si XIV line emission is pro- 
duced, and the Q ratio is smaller. This shows that the 
way the thermal instability is treated in the models (in 
this case, whether we pick the evaporation or condensa- 
tion solutions), has a dramatic effect on all the line ratios 
which are sensitive to the ionization distribution. 

The Hea line triplets can be used as density diagnos- 
tics, but they may be affected by photoexcitation by the 
UV field from the accretion disk (Gabriel & Jordan 1969; 
Blumcnthal, Drake, & Tucker 1972). The forbidden line / 
of O VII at 22.097 A is suppressed due to collisional de- 
population at high density, since n e > 10 14 cm -3 (see the 
density profiles in Fig. 14 and the +7 relative abundance 
distribution in Fig. 15). The O VII intercombination i to 
resonance r line ratio is g = (f + i)/r ~ i/r > 4, indicating 
a purely photoionized plasma (Porquet & Dubau 2000). 
The O VII and N VI Hea line ratios are included in the 
model at their high-density limit, while the Ilea line ratios 
of other ions such as Si XIII have not been modeled yet, 
since the line ratios will start to be a function of position in 
the atmosphere, adding complexity. The depopulation of 
the forbidden line in many He-like ions was attributed to 
resonant photoexcitation in Hercules X-l (Jimenez-Garate 
et al. 2002). The intense UV fields in LMXBs imply the 
same effect will operate (Liedahl et al. 1992). In the con- 
text of the present model, the region where He- like ions are 
abundant is very close (less than 10 7 cm away) to the pho- 
tosphcric surface, such that the UV energy density is as 
high as in the photospheric surface. This implies that the 
He-like diagnostics will be degenerate to high density and 
UV field effects (cf., Mauche, Liedahl, & Fournier 2001). 

The optical depth of the atmosphere may be probed 
by comparing the observed spectra to the model, which 
assumes the lines are optically thin. In particular, the 
O VII r line may differ from the modeled value due to res- 
onant scattering of continuum photons. Whether r gets 
enhanced or absorbed depends on geometry and the rela- 
tive placement of emissivity and opacity. The i and / lines 
should be optically thin, so optical depth can modify the 
g ratio substantially from g ~ 4.2, the value expected for 
a photoionization-dominated, optically thin gas. The Lya 
line in hydrogenic ions also has a large scattering cross 
section. Thus, the Hea r and/or Lya lines are good indi- 



cators of optical depth. 

The RRC can be used for temperature diagnostics, and 
also for probing the behavior of the thermal instability. 
The local RRC width oc T (Liedahl & Paercls 1996). 
The O VII RRC broadening is Doppler-dominated, re- 
sembling the RR lines, for both the evaporating and con- 
densing cases. The O VIII RRC shape varies noticeably 
from the evaporating to the condensing condition. In the 
evaporating case, the O VIII RRC has two temperature 
components, one with a FWHM ~ 2 A that is pro- 
duced at T ~ 10 6 K, and another narrow component with 
T <~ 10 5 K. The two components are distinguishable be- 
cause intermediate temperatures are thermally unstable. 
Thus, the RRC profile in this case provides evidence for 
the existence of this thermal instability. In the condensing 
case, the RRC has a single, broad temperature component 
at T ~ 10 6 K, since the disk atmosphere temperature sud- 
denly drops from the latter value to one where the X-ray 
emission is negligible. The broadening of this RRC is also 
peculiar, since RRC are usually narrow for all photoionized 
gases, given that the kinetic energy of the gas particles is 
generally much smaller than the recombining photon en- 
ergy (Liedahl & Paerels 1996). 

7.1. Luminosity Dependence 

The disk model was run with a lower central luminosity 
L = O.lLEdd, for comparison to the LEdd case explored 
in the previous sections, to investigate the structural and 
spectral changes of the disk. 

The atmospheric radiative recombination luminosity 
was ~ 20 times lower than the L = LEdd case, indicating 
a nearly linear dependence of the disk luminosity to the 
central luminosity. Otherwise, the low-luminosity spec- 
trum shown in Figure 24 shows much resemblance to its 
^Edd counterpart. The dependence of the recombination 
luminosity can be approximated by 

L u ^i cx L x — (19) 

47T 

where fi is the solid angle subtended by the disk atmo- 
sphere. This relationship is then modified by the changing 
density, opacity and thickness of the atmosphere. To ex- 
plain the recombination luminosity behavior, we will first 
describe the calculated atmosphere structure. 

The photospheric and atmospheric boundaries were fit- 
ted by power laws. The fit parameters defined in section 
6 were C ph ot = (1.8 ± 0.3) x 10~ 3 cm, n phot = 1.14 ± 0.01, 
C atm = (1.2 ± 0.2) x 10~ 3 cm, and n atm = 1.18 ± 0.01. 
The fit errors are shown, while systematics are expected 
to follow the same trends as in section 6. The solid angle 
subtended by the entire disk Q/Att ~ 0.1, while for the 
Eddington Luminosity case Si/An ~ 0.2, as can be seen 
by comparing Figure 9 with Figure 25. This implies little 
variation of the disk shape with luminosity. With a factor 
of ten reduction in luminosity, the radiative energy inci- 
dent on the disk is 20 times smaller, which coincides with 
the observed reduction in the recombination emission. 

The accretion disk structure for O.lLsdd, shown in Fig- 
ures 10, 11, 13, and 14, yields a density ~ 4.5 times smaller 
than the LEdd case. Aside from the density change, the 
ionization structure remains quite similar to the LEdd case 
(see Fig. 17 [a] and 17[b]). Naively, a factor of 10 decrease 
in the density would be expected to keep £ constant. How- 
ever, since the atmospheric volume shows little change, the 
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observed density change implies a factor of <~ 20 decrease 
in the ion emission measure, verifying the consistency of 
the density with the modeled recombination flux. 

We reconcile £ with the larger-than-expected density of 
the atmosphere by accounting for a decrease in the atmo- 
spheric opacity. Considering £ to be constant, such that 



X 



(20) 



n n 
oc n, implies that a decrease in density de- 



where r(n) 

creases the opacity, which increases the local flux in 
the atmosphere. Thus, the overdensity of the atmo- 
sphere is explained by a factor of ~ 2 decrease in e~ T , 
which was verified in the models. For the L = O.lLEdd 
case, the atmosphere transmits 34 % of the incident 
flux directly for a disk annulus with r = 10 11 cm, i.e. 
/ F v (z p hot)di'/ J F v (z COI )du ~ 0.34, compared to 12 % in 
the L = LEdd case (see section 6). The atmospheric opac- 
ity is a function of radius, but the observed spectra are 
weighted towards the largest radii. 

The argument above holds for plane parallel atmo- 
spheres that are photoionized by radiation incident at a 
small grazing angle 6 <C 1, such that r ~ 1 for a grazing 
ray, but rCl for any other ray. Therefore, the recombi- 
nation luminosity will be proportional to the total emission 
measure for most viewing angles. 

As more of the X-ray continuum is transmitted onto 
the disk photosphere, the heating of the optically thick 
disk increases. This effect is taken into account to derive 
a self-consistent midplane disk temperature, to 10 %. The 
input grazing angle and the grazing angle calculated from 
z a tm are consistent at the 25 % level. 

7.2. A weak coupling between the disk and its atmosphere 

There is a negligible change in the atmospheric structure 
by varying a from 0.1 to 1 in the disk model. As explained 
below, this can be understood as a decoupling which ex- 
ists between the optically thick disk and the photoionized 
atmosphere. 

The viscosity parameter a has no effect on the disk tem- 
perature, but it does on the density, since p oc a~ l (SS73). 
Assuming vertical isothcrmality in the optically thick disk 
(below the photosphere), its pressure can be obtained from 
equation (3) and is given by: 

P(z) = P e- z2/2Z * (21) 
with P Q oc oT x . The photoionized atmosphere is placed 
on top of this disk, and satisfies two boundary conditions: 
matching temperature and pressure at the photospheric 
boundary. Assume the base of the atmosphere is at some 
pressure P p hot 7 and let a change in viscosity from a to 
a' produce a change in the subphotospheric disk pressure 
from P Q to Pq. To match the boundary conditions, the 
atmosphere has to shift in height, such that 

Pphot = P e- z2/2 ^ = P' e-V?l™l (22) 
which implies that the height shift in the atmosphere, 
Az = z' — z, is given by: 



(23) 



where the approximation z 2 >• Zp was used, which is valid 
since z p hot > 3Zp in the models for a < 1. In such a case, 




it follows that Az/z <C 1, so that very large changes in 
viscosity only produce minute shifts on the atmospheric 
height and negligible effects on the atmospheric emission 
(as already pointed out by Nayakshin, Kazanas & Kallman 
(2000)). 

There are, however, other situations where the role of 
dissipation in the atmosphere has to be reassessed: 1) if 
the dissipated energy in the disk is no longer negligible 
compared to the exterior illumination energy, as might be 
the case for disks around black hole candidates, and if en- 
ergy is transported to the atmosphere via magnetic flares, 
for example, and 2) if there is negligible dissipation, but it 
is sufficient to enhance mixing and, therefore, change the 
atmospheric structure. 

Rozahska, Dumont, Czerny, & Collin (2002) noted that 
the structure of the optically-thick disk has non-negligible 
effects in the corona when the disk scale height is compa- 
rable to or larger than the coronal scale height (in our no- 
tation, when Zp > z cor ). This does not apply to our case, 
because for a neutron star LMXB disk with r > 10 8 5 cm, 
we get z cor S> Zp. Rozahska, Dumont, Czerny, & Collin 
(2002) assert that the disk and corona are coupled in BHC 
for r ~ 10 Rq, where Rq is the Schwarzschild radius; while 
for AGN, the situation depends on the accretion rate. By 
this measure, the disk-corona coupling may also be signif- 
icant in the outer radii of BHC disks, if the reduction of 
the central illumination flux shrinks the size of the corona. 

7.3. Emission line profiles 

Synthetic profiles were produced for all the emission 
lines. For simplicity, the calculation assumes that line 
scattering is negligible, which may not be a valid assump- 
tion for resonance lines. Here we select the profiles of the 
brightest lines that are not contaminated by other ions. 
The line profiles are also calculated for a disk within the 
10 8,5 < r < 10 10 cm radius range, for comparison with 
the 10 8 5 < r < 10 11 cm radius range shown above. The 
emission line profiles have more broadly separated peaks 
for smaller disks, as shown in Figure 26. 

No variation in line broadening is obtained as a function 
of charge state. The N VI Lya line has the same velocity 
profile as the corresponding lines in O VIII, Ne X, Si XIV, 
and Fe XXVI. We attribute this to the vertical stratifica- 
tion of the atmosphere, which allows the full range of ion- 
ization parameters detectable in the X-ray band to exist 
in every annulus. Future two-dimensional models taking 
into account radiation transfer in the radial direction may 
exhibit a trend for the line widths, since the radial optical 
depth is not negligible. 

8. DISCUSSION 

8.1. Neutron star LMXB spectra observed with Chandra 
and XMM-Newton 

The LMXB spectra observed with the Chandra High 
Energy Transmission Grating (HETG) and the XMM- 
Newton RGS already provide stringent tests for the mod- 
els in this article. The high resolution spectra of LMXBs 
are generally dominated by continuum emission, which 
is sometimes punctuated by emission or absorption lines. 
Only a fraction of the sample of observed neutron star 
LMXBs show prominent lines. The accreting pulsar 
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4U1626-67 shows double-peaked and broad emission lines 
(Schulz et al. 2001). The eclipsing dipper EXO0748-67 
has broad emission lines (Cottam et al. 2001a). The ac- 
cretion disk corona source 4U1822-37 has narrow emission 
lines (Cottam et al. 2001b). The dipper source 4U1624- 
49 shows narrow absorption lines (Parmar, Oosterbroek, 
Boirin, & Lumb 2002). Her X-l, an intermediate-mass 
X-ray binary with a precessing accretion disk, has narrow 
emission lines during its low- and intermediate-flux states 
(Jimcnez-Garate et al. 2002). The line emission spectra 
are dominated by H-like and He-like ions. P Cygni pro- 
files are observed in Circinus X-l (Brandt & Schulz 2000), 
a unique LMXB with high- velocity outflows. Most LMXBs 
do not exhibit discrete spectral features aside from inter- 
stellar absorption, such as X0614+09 (Paerels et al. 2001). 
The emerging pattern implies that at least three classes of 
neutron star LMXBs exist which produce detectable X-ray 
emission lines: 1) high inclination LMXBs with i > 70°, 2) 
accreting X-ray pulsars, and 3) LMXBs with high velocity 
winds. 

Our models are proving to be of great relevance to the 
interpretation of the spectra of LMXBs with high incli- 
nation and LMXBs with an X-ray pulsar. Most impor- 
tantly, all LMXB spectra validate the basic assumption 
in our model: the plasma is photoionized. Evidence for 
other heating mechanisms, such as shocks, is not observed 
in LMXB spectra. Shocks are predicted by the Miller & 
Stone (2000) MHD disk models, but it is not clear whether 
X-ray emission from such shocks would be observable. The 
observed signatures of a plasma heated primarily by pho- 
toionization are the RRC, the peculiar Hea line ratios, 
and the weakness of Fe L line emission relative to that 
of low-Z and mid-Z elements (Liedahl 1999). The spectra 
have prominent line emission from H-like and He-like ions. 
These properties are shared by all the spectra shown in 
this article. Furthermore, the line velocity broadening ob- 
served in two LMXBs (4U1626-67 and EXO0748-67) pro- 
vides kinematic evidence for accretion disk atmospheric 
and coronal emission. 

By contrast, our model is inadequate for the interpre- 
tation of the spectra of LMXBs with high- velocity winds. 
These spectra show a photoionized gas with significant op- 
tical depth. Circinus X-l is rare because of the extreme gas 
dynamics and the sizable line optical depths which are evi- 
dent in its spectrum. The density n e ~ 10 13 cm~ 3 deduced 
from the observed P Cygni profiles in Circinus X-l (Brandt 
& Schulz 2000) is consistent with our model calculations, 
but the wind dynamics rules out our assumption of hydro- 
static equilibrium. The modeling of an LMXB disk wind 
spectrum requires a wind acceleration mechanism and a 
revised disk structure (Proga & Kallman 2002). 

Initially, our model can be used as a tool for the identi- 
fication of the discrete X-ray spectral signatures from the 
accretion disk atmosphere and corona. Our model pro- 
vides a quantitative expectation of the X-ray line fluxes 
produced by the entire disk. Our physical calculation of 
the line profiles (and therefore the line emissivity as a func- 
tion of r), can be used to measure the maximum disk ra- 
dius and the radial ionization distribution. Furthermore, 
we have modeled the density at which each ion is pro- 
duced, and this is testable with plasma diagnostics. The 
ionization distribution (measured by the flux of the emis- 



sion from the high-ionization species relative to that of 
the low-ionization species) is a fingerprint of the disk at- 
mosphere and corona. The temperature structure of the 
disk atmosphere can also be probed with the temperatures 
measured with the RRC of various ions. A full investiga- 
tion of the models as compared with spectral data will be 
performed in a future paper. 

8.2. Limitations of the model 

The disk structure can be improved by relaxing the 
assumptions made in the radiation transfer calculations. 
The spectral model calls for the inclusion of additional 
lines. The observed spectra will also allow us to investi- 
gate additional physics in the disk which may be missing 
in our current model. 

A 2-D or 3-D transfer calculation is needed to improve 
the coronal structure model. To simplify the radiation 
transfer calculations, we split the disk into a set of nested 
cylindrical shells, and we use 1-D transfer to calculate the 
structure of each shell. By defining z atm (r) and iterat- 
ing on 8(r) in each shell until self-consistency is obtained 
(see eq. [8] and Fig. 5), we produce a pseudo-2-D trans- 
fer calculation. However, the 1-D transfer approximation 
starts to break down at the largest radii, since z cor ~ r at 
r > 10 10 cm (see Fig. 10). 

The UV emission is included in the structure model but 
remains to be added to the high resolution spectral model. 
However, we do not expect substantial differences between 
the UV spectrum obtained with our disk model and the 
results by Raymond (1993). The optical depth of UV lines 
such as C IV is ^ 100 in the latter model. The UV and op- 
tical lines originate just above the photosphere (z p hot), at 
densities of 10 13 < n e < 10 cm~ 3 . The structural differ- 
ence between our models and those by Raymond (1993) 
occurs at the X-ray emitting atmosphere and corona. In 
contrast, the z p hot(f) m our models agrees well with those 
by Raymond (1993). 

Fluorescence and resonant scattering will need to be 
added in the spectral model. The fluorescence line flux 
can be of the same order as the recombination line flux. 
The r-dependence of the Fe K fluorescence flux should be 
distinct from the r-dependence of the recombination emis- 
sion. A 6.4 keV Fe K emission line is produced by M-shcll 
charge states of Fe localized at the base of the atmosphere 
(see Fig. 17). The Fe K fluorescence flux will scale with the 
hard X-ray transmittance of the atmosphere and corona. 

The structure model indicates that the optical depth of 
resonance lines is large. However, a realistic treatment of 
resonant scattering in our spectral models is complicated 
by the velocity shear within the disk. Resonant scattering 
and Fe K fluorescence emission were only included in the 
low-resolution 1-D transfer calculations to obtain the disk 
structure. The propagation of resonance line photons is 
highly anisotropic, and it depends on the viewing angle 
because of the Keplerian velocity shear and the geomet- 
rical thickness of the atmosphere. LMXBs with strong 
emission lines do not exhibit detectable features from res- 
onant scattering of continuum photons, with the exception 
of Circinus X-l (section 8.1). 

Magnetic fields, which arc not included in our model, 
may affect the structure of some regions of the disk (see 
section 8.3). 



14 



Jimcncz-Garate et al. 



8.3. A strongly magnetized corona? 

A strong magnetic field may affect the coronal structure 
and the X-ray spectrum. The corona under considera- 
tion is located at the outer radii of a centrally illuminated 
disk. We believe B fields play a secondary role in this 
type of corona, because the energy budget of the corona 
is dominated by X-ray irradiation. In detail, the role of B 
fields cannot be discounted, because phenomena such as 
magnetic flare heating may dominate over photoionization 
within localized regions of the disk. Recent MHD models 
predict B fields larger than the virial value in the disk 
corona. However, the applicability of these MHD models 
to the illuminated corona is dependent on the effects of 
radiative heating and magnetic reconnection. 

MHD models of radiationless accretion disks show that 
above a few scale heights, the magnetic pressure is larger 
than P gas (Miller & Stone 2000). The gas dynamics in 
the disc is dominated by the MRI. Miller & Stone (2000) 
found that the MRI produces B fields which buoyantly 
rise to the atmosphere and corona. In their model, 25% 
of the magnetic energy generated by the MRI rises to the 
corona, representing 60% of the local heating, but < 4% 
of the dissipative heating in the disk. A 3D MHD disk 
model by Machida, Hayashi, & Matsumoto (2000), with 
an initial toroidal configuration, shows that the strong B 
fields in the corona are confined in filaments, with a filling 
factor of a few percent. Another 3D MHD model by Haw- 
ley, Balbus, & Stone (2001) confirms the presence of large 
B fields in the corona. 

By contrast, when X-ray illumination is present, the B 
field plays a relatively minor role in the overall energetics 
of the disk corona. This is true at least in a spatially- 
averaged sense. The magnetic energy produced by the 
MRI can be no larger than the energy dissipated in the 
disk. The maximum energy available for the B field scales 
as r~ 3 , and it is given by the first term on the right side 
of equation (1). Assuming all of the accretion energy is 
contained in the B field, equation (2) indicates that the 
illumination energy is larger than the magnetic energy for 
r > 10 10 cm (per unit disk area). Since MHD disk mod- 
els typically assume an isothermal or adiabatic disk, they 
may not apply to an extended hot corona dominated by 
photoionization. 

Magnetic reconnection has the effect of decreasing the 
magnitude of the disk-coronal B field. Coronal flares 
resulting from magnetic reconnection convert magnetic 
energy into kinetic energy through particle acceleration. 
Most of the X-rays in black hole accretion disks may be 
produced by magnetic flares (di Matteo 1998). In the Liu, 
Mineshige, & Shibata (2002) model, reconnection events 
in disk flares reduce the B field by a factor of ~ 30. Thus, 
the buoyant B field may be dissipated in a flaring region. 

9. CONCLUSIONS 

We have calculated the hydrostatic structure of a pho- 
toionized accretion disk atmosphere which is in thermal 
equilibrium and ionization balance. We also determined 
the atmosphere's thermal stability and its observable high 
resolution X-ray recombination emission spectrum. 

• A feedback mechanism between illumination and 
atmospheric structure enlarges the atmosphere. 



The disk atmosphere is orders of magnitude less 
dense than the disk midplane. The atmosphere 
extends for a few tens of disk pressure scale heights 
(if the pressure scale height is calculated using the 
disk midplane temperature). Illumination heats 
and expands the disk atmosphere, increasing the 
number of absorbed photons in the atmosphere 
and heating it further, producing further expansion 
of the atmosphere, and so on. The expansion 
stops because the atmosphere becomes optically 
thin, cooling and contracting. The inclusion of 
the feedback mechanism increases the size and the 
line emission flux of the atmosphere by an order 
of magnitude. The atmospheric thickness is much 
larger than the standard a-disk model thickness, 
and it is consistent with the <~ 12° subtended 
semi-angle deduced from optical modulations in 
LMXBs. The disk atmosphere thickness also 
explains the under- abundance of eclipsing LMXBs. 

• The disk atmosphere subtends a large solid 
angle 0.07 < tt/Air < 0.2. If the inclination is 
i > 80°, the disk photosphere (which subtends 
0.04 < VL/Ak S 0-08) may shield the neutron 
star flux, producing an ADC source with partial 
eclipses or without eclipses altogether. The disk 
f2 depends weakly on the neutron star luminosity, 
but f2 scales linearly with disk radius. The disk 
recombination luminosity scales linearly with f2. 

• The atmospheric structure is independent of the 
viscosity parameter a. The viscosity changes the 
density in the optically-thick part of the disk, 
producing a small shift in atmospheric height, but 
this has no effect on the X-ray spectrum. 

• The X-ray spectra are dominated by lines from 
H-like and He-like ions of abundant elements from 
C to Fe, as well as RRC and weak Fe L lines. The 
line ratios are a sensitive probe of the atmospheric 
and coronal structure. 

• Clear spectral signatures of photoionization are 
present, as well as temperature, density, and 
radiation field diagnostics. An intercombination 
to resonance line ratio of ~ 4 is modeled 

for low-Z He-like ion line triplets. RRC are 
unequivocal signposts of photoionization. The 
density diagnostics from Hea lines of low-Z and 
intcrmcdiate-Z elements are degenerate with the 
presence of an intense UV radiation field from the 
disk itself, so the R ratio may not give conclusive 
signatures of high-density in LMXBs. Much of the 
disk atmosphere is close to the photosphere, such 
that the dilution factor of the UV field is small. 
The Hea density diagnostics could operate at the 
densities predicted by the disk atmosphere model 
in LMXBs, of 10 13 < n e < 10 15 cnr 3 . 

• The line fluxes are nearly proportional to the 
X-ray continuum luminosity. The disk line fluxes 
decreased by a factor of 20 when the system 
luminosity was decreased by a factor of 10 (to 

L = 10 37 - 3 erg s 1 ). The atmospheric density 
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was reduced by a factor of ~ 5, its optical 
depth was reduced, and the atmosphere was ~ 2 
times less extended than in the high-luminosity 
(L = 10 38 - 3 erg s" 1 ) case. 

• The line equivalent widths depend strongly on 
inclination. The relative obscuration of the 
neutron star affects the equivalent width and 
detectability of the disk X-ray emission. The 
modeled disk emission is almost undetectable 
when the neutron star continuum is also in the 
line of sight. As such, high inclination systems, 
or systems with dips or ADC, are more likely to 
show X-ray lines due to enhanced contrast. This 
expected trend has been largely confirmed by 
Chandra and XMM-Newton observations. We have 
demonstrated that for a highly absorbed neutron 
star continuum in our fiducial system, the disk 
X-ray lines are detectable with both the Chandra 
and XMM-Newton grating spectrometers. 

• Double-peaked X-ray lines can be detected for 

r = 10 10 cm disks, but larger r = 10 11 cm disks 
may appear blended in a single peak in the grating 
spectra. The emission line region spans several 
orders of magnitude in disk radius. The modeled 
line profiles arc needed to deduce the outer disk 
radius. Line emission from the outer regions of the 
disk dominates. The emission increases with disk 
radius, and the Doppler broadening of the lines 
decreases for larger r. The wings of the broadest 
lines are lost in the continuum, decreasing their 
apparent equivalent width. 

• The resonance line optical depths can be measured. 
If the r line in He-like ions has a value which 
differs from the model calculations, it may be due 
to resonant scattering of continuum photons. The 
line ratios in the Lyman series can also work as 



optical depth diagnostics. We have not included 
these effects on the current version of the model, 
but our results for ionic column densities show 
that appreciable line optical depths are present, 
and hence this process will be included in future 
versions of the code. 

• The continuum optical depth of the atmosphere is 
generally small (t <C 1), except for photons which 
propagate nearly parallel to the disk plane. The 
atmosphere is optically thick to X-ray continuum 
photons from the neutron star. However, most of 
the recombination line emission is not appreciably 
affected by the continuum opacity. 

• The spectrum is sensitive to a thermal instability 
present in photoionized gases. By forcing all 
the chosen solutions to be thermally stable, a 
break in the temperature, density, and ionization 
structure is created. Measurably different X-ray 
spectra are obtained depending on the resolution 
of this instability. The shape of RRC profiles, 
which in the models show multiple temperature 
components, and the relative intensity of lines 
such as Mg XI, are useful diagnostics of the stable 
temperature regime. 

The spectra obtained with the Chandra HETG and the 
XMM-Newton RGS show that the plasmas in LMXBs 
are photoionized, as our model assumes. Two LMXBs 
(4U1626-67 and EXO0748-67) show kinematic signatures 
of accretion disk emission (Schulz et al. 2001; Cottam et al. 
2001a). The line fluxes, line profiles, the ionization distri- 
bution, density, and RRC temperatures, provide a wealth 
of diagnostic capability for the identification of accretion 
disk atmospheres and their properties. The spectral com- 
parisons with the data are promising, and they will be 
addressed in future work. 
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Fig. 1. — Schematic of a low-mass X-ray binary (LMXB) system with a neutron star primary. Cross section view. The inclination angle i 
is defined. 
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Fig. 2. — Side view schematic of an illuminated accretion disk and the model geometry, assuming an extended corona above the disk. To 
compute the ionization structure of the disk, the disk geometry (a), can be approximated by a series of gas columns which are illuminated 
from the top (b). 




Fig. 3. — Gas column geometry for each disk annulus. The boundary conditions on the atmospheric structure are shown. 




Fig. 4. — Schematic of the feedback between radiative heating and disk geometry. The heated atmosphere expands and collects more 
radiation, reaching equilibrium at ~ 10 times its initial volume. 
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Fig. 5. — Grazing angles 9(r) for the radiation impinging on the disk, at successive model iterations (A to D). The 8(r) input to the model 
are compared to the 8(r) extracted from the output disk structure. The resulting model D is self-consistent and has a finer grid. 
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Fig. 10. — Modeled temperature (T) vs. normalized disk height (z/r) for various radii. The diagrams can be directly compared to the 
Raymond (1993) result. Our new model atmosphere is 10 X more extended. The chosen stable solutions correspond to a disk under 1) steady 
evaporation (solid line), 2) steady condensation (dotted line), and 3) steady evaporation at low luminosity L = O.lZ/Edd (dashed line). The 
large-scale features of the first two cases are nearly identical. 
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Fig. 11. — Modeled proton density (n p ) vs. disk height (z) for various radii. Evaporating (solid line), condensing (dotted line), and 
low-luminosity (dashed line) disk models, as in Fig. 10. The height is normalized to the local radius. 
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Fig. 12. — Dominant heating and cooling mechanisms versus the vertical height of the atmosphere, for the disk annulus with r = 10 11 cm. 
A Compton-hcated corona and a recombining atmosphere can be discerned. 
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Fig. 13. — Modeled temperature (T) vs. disk atmosphere height above the photosphere (z a tm — Zphot) f° r various radii. Evaporating (solid 
line), condensing (dotted line), and low-luminosity (dashed line) disk models. The height is normalized to the local radius. 
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Fig. 14. — Modeled proton density (n p ) vs. disk atmosphere height above the photosphere (z a tm — z p hot) for various radii. Evaporating 
(solid line), condensing (dotted line), and low-luminosity (dashed line) disk models. The height is normalized to the local radius. 
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Fig. 15. — Vertical distribution of low-Z ions for the evaporating disk model. Low charge states are omitted for clarity. The relative 
ion abundances are plotted against (z a t m — z phot )/r, the vertical height of the atmosphere, for the disk annulus with r = 10 11 cm, which 
dominates line emission. 
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Fig. 16. — Vertical distribution of mid-Z ions. Low charge states are omitted for clarity. The relative ion abundances arc plotted against 



the vertical height of the atmosphere, for the disk annulus with r 
for the evaporating disk model, where they arc strongest. 



10 cm. The thermal instability suppresses the Hc-likc ion lines, even 
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Fig. 17. — Vertical distribution of K and L-shell ions of Fe, plotted for: (a) the condensing disk with L = I/Edd, (b) the evaporating disk with 
L = ^Edd! and (c) the evaporating disk with low-luminosity (O.lZ/Edd)- The relative abundances are plotted vs. the height (z a tm — z p hot)/ r > 
for the disk annulus with r = 10 11 cm. The plateau and break at (z a tm — z p hot)/ r = 0.0001-0.003 correspond to the instability regime. 
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Fig. 18. — Vertical distribution of low-Z ions for the condensing disk model. Low charge states arc omitted for clarity. The relative 



ion abundances are plotted against (z a tn 
dominates line emission. 



Zphot)A, the vertical height of the atmosphere, for the disk annulus with r = 10 cm, which 
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Fig. 19. — Modeled recombination emission from the disk atmosphere (L^ lsk = d 2 E 2 F^ s]i ). Neither the ionizing continuum, interstellar 
absorption, nor Dopplcr broadening are included. 
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Fig. 20. — Model spectrum of an LMXB, with Ai? x = 2 eV bins. The line luminosity from Fig. 19 is Dopplcr shifted by the projected local 
orbital velocity, at an inclination of 75°. For the disk lines, we let jVjj lsk = 10 21 cm~ 2 be constant. For the neutron star continuum, we let 
(a) Nfj = 10 21 cm" 2 , (b) = 5 x 10 22 cm" 2 , and (c) Nfj = 10 24 cm . Each term in eq. (18) is shown, as well as their sum (evaporating 
disk model). 
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Fig. 21. — Spectra for evaporating and condensing accretion disks. Simulated 50 ks observation with XMM-Newton RGS 1. The evaporating 
disk spectrum here corresponds to Fig. 20(b). The continuum emission from the inner (r < 10 8 5 cm) disk is not included, and the obscuration 
of the ionizing continuum from the compact object is assumed. 
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Condensing disk spectrum (as in Fig. 21). Simulated 50 ks observation with Chandra MEG, +1 order. 
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Fig. 23. — Evaporating and condensing disk spectra with occulted neutron star. Simulated 50 ks observations with the Chandra medium 
energy grating MEG, +1 order. The evaporating disk spectrum shown corresponds to the model in Fig. 20(c), with the neutron star 
practically occulted by Nfj = 10 24 cm -2 . 



Evaporating disk model, low luminosity 




15 20 
Wavelength (1) 



35 



Fig. 24. — Model spectrum of an LMXB with L = 0.1L Edd , with A£ x 
is included at an inclination of 75°. The disk emission has Njj = 10 21 cm" 
to the I/Edd case m Fig- 20(b). 
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Fig. 25. — Total solid angle subtended by the disk photosphere (2 p hot)i atmosphere (z a tm), and corona (T = 5 x 10 6 K isotherm). The 
disk structure for the low-luminosity (L = 0.1L Edd ) case is shown, whose atmosphere subtends a solid angle 30 to 40 percent smaller than 
the L = LEdd case on Fig. 9. 




Fig. 26.— From top to bottom and left to right: line profiles for N VI Hea, N VII Lya, O VII Hea, O VIII Lya, Ne X Lya, and Si XIV 
Lya. The bin size is AA = 0.005 A. The solid line corresponds to a disk with maximum radius of r = 10 11 cm, while the dashed line is the 
weaker and broader line profile produced by a disk with r = 10 10 cm. 
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Table 1 

Line fluxes for disk models. Line fluxes in units of 1CT 12 erg cm -2 s _1 , typically over an interval E/AE ~ 100. 
Estimated systematic normalization error limit of ~ 50 %, due to the 1-D transfer calculation. These line 
fluxes do not include any interstellar absorption effects, unlike those from Table 1 in Jimenez-Garate, 

Raymond, Liedahl, & Hailey (2001). 



L — O.lLEdd L — L^dd L — £Edd L = Lftdd 

r = 10 11 cm r — 10 11 cm r = 10 10 cm r = 10 11 cm 



Linc(s) 
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Evaporating 


Evaporating 


Condensing 
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APPENDIX 
RADIATIVE RECOMBINATION EMISSION 

We describe the numerical calculation of the radiative recombination emission, including both recombination lines and 
radiative recombination continua. A similar method for an optically thin gas in the photoionized wind of a High Mass 
X-ray Binary HMXB was described by Sako et al. (1999). 

Consider an infinitesimal volume dV at which a single ionization parameter £, temperature T, electron density n e , and 
elemental abundances A z describe the state of a gas. The £(T) function is found from thermal balance and ionization 
equilibrium, for a given ionizing spectrum F v (section 3.3). In the radiative recombination process 

Z +{1+1) + e" -> Z+" + hv RRC (Al) 

an electron recombines with an ion with net charge +(i + 1), assumed in its ground state, producing a new ion with net 
charge +i, which might be excited. The radiative recombination continuum photon has energy 

£ x = hv RR c = x + KE e , (A2) 

where \ is the ionization energy of ion Z +l * , and KE e is the initial kinetic energy of the electron, assumed to be on a 
Maxwell distribution with temperature T. The radiative recombination rate to Z +l (where Z +l can be in any quantum 
state) in units of s _1 is 

Tflfl = n c n z ^ +1 a RR dV (A3) 

where n z ,i+i is the Z + ( l+1 ' ) ion number density, equation (A3) defines ozrr, the total radiative recombination rate coefficient 
in units of cm 3 s _1 . Note cxrr depends on Z, i, and T. 

Recombination Lines 

After recombination, a fraction r\ u ^i of the Z +l ions produce a radiative cascade photon by an electronic transition 
from upper level u to lower level I. The line luminosity of photons from this transition in units of erg s _1 is 

dL u ^i = n e n z , i+1 E u ^ir] u ^iaRRdV (A4) 

where E u ^i is the transition energy in ergs. To create a synthetic spectrum, the line luminosities dL u ^i at energies E u ^i 
in the X-ray band are added for all the levels u,l, and all the ions Z + ( l+1 ^ which are abundant in the gas. Notice the 
recombination emission of the Z +l ion depends on the number density of the Z + ( l+1 ^ ion. 

For computational purposes, various quantities from equation (A4) are defined. The specific line power 

S u ^i = Vu^iOiRR, (A5) 

in units of cm 3 s _1 , is the photon emission rate per ion, per unit electron density. The Z +( - t+1 ' ion was assumed 

to be in its ground state before recombining into Z +l . The population fraction f u of each level of the Z +% ion is computed 
explicitly, and S u ^i is obtained by equating the matrix of the photon emission rates per ion, 

n c S u ^i = f u A u ^i (A6) 

where A u ^i is the rate of spontaneous decay for Z +l . After solving S u -,i for a grid of temperatures, typically in the 10 
to 80 cV range, it is fit to a power law 

Su^i = C U ^T-^' (A7) 
where the exponent 7 U _>; is typically 0.6-0.8. The number density of Z + ( l+v> is calculated with 

n z ,i +1 = n H A z f Zii+1 (A8) 

where njj is the proton density, A z is the fractional abundance of element Z relative to H, and f z .i+i is the fractional 
abundance of the Z + ( l+1 ^ ion relative to all the Z ions. The differential emission measure for is defined as 

d(EM Zti+1 ) = n c n Ztl+1 dV, (A9) 

in units of cm~ 3 . The line luminosity in equation (A4) can therefore be re- written as 

dL u ^ = E u ^ l S u ^id{EM Z!i+1 ). (A10) 

If the emission measure is defined as d(EM) = n^dV, then dL u ^i = P u ^id(EM), where P u ^i is defined as the line 
power, with units of erg cm 3 s _1 . The emission measure is useful for calculating the luminosity of an optically thin gas. 
Since the accretion disk atmosphere does have some optical depth, d(EM) and d(EM z ^ + i) will only be used to track the 
regions where the emission originates. The radiative recombination line list includes transitions from levels with principal 
quantum number n < 4 or 5 typically, although in some cases levels of up to n = 7 are included. 
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Radiative Recombination Continuum 

To calculate the shape and luminosity of the RRC, a Maxwell thermal distribution, the photoionization cross sections, 
and the Milne relation were used. The monochromatic version of the RR coefficient in equation (A3), for electrons with 
velocities between v and v + dv is 

olrr, v = o- RR , v vf v dv, (All) 

where orr, v is the RR cross section of ion Z + ( l+1 \ and the number of electrons in that velocity range is f v dv, which is 
assumed to be given by the Maxwellian distribution 



1/2 , \ 3/2 



f * = {5) \w ^- mv2/2kT ' ( A12 ) 

where m is the electron mass. Thus, the monochromatic RRC emissivity of Z +l for thermal electrons is 

dv 

j v =n e n Zii+ iE x a RRtU vf v -p=-, (A13) 

in units of erg cm~ 3 s _1 erg -1 . Because radiative recombination is the inverse process of photoionization, a relationship 
between their cross sections is derived by equating their transition rates obtained from Fermi's golden rule (Salzmann 
1988). Detailed balance yields a cross section ratio proportional to the ratio of the density of final states for each reaction. 
For recombination and photoionization, this is the Milne relation 

9i ( E x 



0-RR,v = 0-PEm, (A14) 

9i+i \mcvj 

where ope,v is the photoionization cross section for the valence electron of Z +l , and gi,g i+ i are the statistical weights of 
the energy levels of ions Z +l and Z + ( l+v> , respectively. Note g = 2 J + 1, for total angular momentum quantum number 
J. From equation (A2), and equations (A12)-(A14), one can derive the RRC emissivity 

/ 9 \ 1/2 „ / P 2 \ 3/2 

>=(-] n c n z , +1 ^ca PE J-%-) e~^-^ kT (A15) 



/ 9i+i \mc kl 

which is in the same units as equation (A13). The ground state photoionization cross sections are taken from Saloman 
Hubble, & Scofield (1988), and the code based upon that paper is used to calculate the cross sections from excited levels 
such that 

Rv r 3 / ^ x q 
opE.v = lO^V-^exp 
X 



(A16) 



L 9 =o 

in units of cm 2 , where the four-element a q vector and \' are fitting parameters, and Ry = 13.6 eV. Note x' ~ X- The 
constant n' is a function of various occupancy numbers and statistical weights. 

CONTINUUM OPACITY 

In the disk atmosphere, the recombination emission is partially absorbed by the ionized gas above it. Each ionization 
zone in the gas column in Figure 3 is denoted by an index j = 1...N, starting from the top zone. If a recombination 
emission net flux F v j is produced in each zone j of height hj , then the total flux for the column is 

N 



: ^ h r , 



COS I 



(Bl) 



F v = F ",3 ex P 

i=i 

where i is the inclination angle of the observer in reference to the disk midplane normal, and K„ im is the continuum opacity 
of the mth zone, 

Kv,m = 0-Tn e ,m + ^ O v,z.kn z ^,m , (B2) 
z,k 

where n z ^, m is the number density of each ion Z +k in the mth zone, n e , TO is the electron density, ax is the Thomson 
cross section, and <r„ )Z) k is the photo-electric absorption cross section of ion Z +k , given by 

z — k 

0~v,z,k — Q~PE,v,z.k,e (B3) 

e=l 

where the photoionization cross section o-pE,u,z,k,e for each electron e in the ion Z +k is given by equation (A16), and 
the cross sections for all Z — k electrons were added. This contrasts with the case of recombination, in which only the 
vpE,v for the valence electron was needed. The model atmosphere is optically thin; i.e., the continuum optical depth 

t v = hjK V j <c i. 

The flaring geometry of the disk atmosphere (see section 6) is not taken into account in equation (Bl). Disk flaring 
will result in opacities that are larger than those in equation (Bl) at inclinations i > 75-80°, since the disk atmosphere 
subtends an angle of arctan (z atm /r) — 10-15°. 



